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Abstract: Starting from the strong-coupling SU(2) Wilson action in D = 3 dimen- 
sions, we derive an effective, semi-local action on a lattice of spacing L times the 
spacing of the original lattice. It is shown that beyond the adjoint color-screening 
distance, i.e. for L > 5, thin center vortices are stable saddlepoints of the corre- 
sponding effective action. Since the entropy of these stable objects exceeds their 
energy, center vortices percolate throughout the lattice, and confine color charge in 
half-integer representations of the SU(2) gauge group. This result contradicts the 
folklore that confinement in strong-coupling lattice gauge theory, for D > 2 dimen- 
sions, is simply due to plaquette disorder, as is the case in D = 2 dimensions. It also 
demonstrates explicitly how the emergence and stability of center vortices is related 
to the existence of color screening by gluon fields. 
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Instantons 



Quark confinement is commonly attributed to the influence of some special class 
of gauge field configurations, which dominate the QCD vacuum at large scales. Be- 
cause of their high probability, these dominant configurations would most naturally 
correspond to the saddlepoints of an infrared effective action, derived at large scales 
by integrating over high-frequency modes. In strong-coupling lattice gauge theory 
there are methods available which enable us to compute the QCD spectrum and 
string tension analytically, and the same methods could also be applied to extract 
a long-range effective action. An interesting question, then, is what type of saddle- 
point configurations are actually found at strong lattice couplings; it is likely that 
the answer would also shed some light on QCD in the continuum limit. 

The classical Euclidean action of pure SU(N) gauge theory is stationary, or nearly 
so, at multi-instanton configurations. In quantized lattice gauge theory, however, 
we can imagine performing renormalization-group (RG) transformations so as to 
obtain an effective action at some scale R. For scales R well below the confinement 
scale, the main effect of the RG transformations will simply be the running of the 
lattice coupling constant. At larger scales, however, so-called irrelevant operators 
can become important in the effective action. As a consequence, at these larger 
scales, the effective theory may have non-trivial saddlepoints which are something 
other than instantons. 

There are good reasons to believe that at sufficiently large scales, these non-trivial 
saddlepoints are center vortices. On the theoretical side, we note that the asymptotic 
string tension between static color charges in SU(N) gauge theory depends only on 
the iV-ality of the color charge representation. Although this fact is deduced rather 
trivially from the possibility of color-screening by gluon fields, it has some profound 
implications for the infrared structure of the QCD vacuum. Consider, for example, 
Wilson loops Wj{C) in SU(2) gauge theory, where j = 0, |, 1, |, ... labels the group 
representation. Wilson loop expectation values can be viewed as a probe of vacuum 
fluctuations in the absence of external sources (think of evaluating a spacelike loop 
in the Hamiltonian formulation), and large Wilson loops are presumed to become 
"disordered," i.e. have an area-law falloff, due to averaging over certain types of large- 
scale fluctuations which dominate the vacuum state. Whatever the nature of these 
confining fluctuations, they must have the highly non-trivial property of disordering 
only the j = half-integer loops, but not the j = integer loops. Center vortices are 
the only configurations we know of that have this property. Dual-superconductor 
models, in which all multiples of abelian electric charge (identified in an abelian- 
projection gauge) are confined by the dual Meissner effect, do not seem satisfactory. 
In these models, the potential between charged objects is roughly proportional to the 
electric charge. But this charge dependence cannot be correct, since gluons carrying 
two units of electric charge are available to screen multiply charged sources, and 
numerical simulations indicate that only odd multiples of the abelian charge (non- 
zero iV-ality) are confined, while even multiples of abelian charge (zero iV-ality) are 
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screened [1]. Center vortices seem to be the natural way of accounting, in terms 
of dominant field configurations, for this very fundamental distinction between zero 
and non-zero iV-ality, which is evident even in the abelian projection. 

On the numerical side, there is now abundant evidence in favor of the vortex 
theory of confinement [2-13], which we will not attempt to review here. The present 
situation can just be summarized as follows: There exists a method (known as "center 
projection") for locating center vortices on thermalized lattices; the rationale unde- 
lying this method is explained in ref. [14]. By locating the vortices, their effects on 
gauge-invariant observables such as Wilson loops, Polyakov lines, topological charge, 
etc., can be studied in detail. The numerical evidence indicates that fluctuations 
in vortex linking number are the origin of the asymptotic string tension of Wilson 
loops. The free energy of a vortex, inserted into a finite lattice via twisted boundary 
conditions, has also been computed, and has been shown to fall off exponentially 
with the lattice size at just the rate predicted by the vortex theory [12, 13]. 

Vortices presumably have a finite thickness comparable to the adjoint string- 
breaking length, at about 1.25 fm [15], where the crossover from Casimir scaling to 
iV-ality confinement occurs (cf. the discussion in ref. [16]). Independent estimates of 
the vortex thickness are based on measurements of "vortex-limited" Wilson loops in 
ref. [4], and on the vortex free energy in finite volumes [12]. Both of these estimates 
put the vortex thickness at a little over one fermi. Beyond this scale, the presence of 
vortex sheets in the vacuum should be very evident. A reasonable conjecture is that 
if the appropriate effective action could be determined at this scale, it would be found 
to have stable saddlepoints corresponding to vortex configurations, which percolate 
through the lattice according to the usual energy-entropy arguments. Unfortunately, 
the calculation of long-range effective actions is very difficult even numerically, via 
Monte Carlo RG methods, and at such large scales the problem is quite intractable 
by perturbative (e.g. one-loop) methods. 1 

We therefore turn our attention, in this article, to strong-coupling lattice gauge 
theory, where analytic methods can be brought to bear at arbitrarily large scales. In 
the strongly coupled theory in D > 2 dimensions, we have a linear static potential for 
all color charge representations up to a screening scale of about four lattice spacings 
(for the SU(2) Wilson action). Beyond that scale, the string tension depends only 
on the iV-ality of the representation, just as in the continuum theory. Thus, if our 
conjecture is correct, and if the iV-ality dependence implies a vortex mechanism, the 
effective strong-coupling action at a scale beyond four lattice spacings should have 
saddlepoints which are stable center vortices. 

1 On the other hand, there do exist some intriguing results at one- loop that should be noted. 
Diakonov [17] has computed a one- loop effective potential for magnetic flux tubes, and his result 
indicates that the potential is minimized for magnetic flux in the center of the gauge group. There 
is also the old, but still provocative, Copenhagen vacuum picture, which is again based on one-loop 
considerations [18]. 
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Figure 1: Decomposition of a large area bounded by planar loop C into many sub-areas 
bounded by planar loops {Cj} 



There is, however, some folklore to the effect that confinement in strong coupling 
lattice gauge theory, in D > 2 dimensions, is essentially the same as in D = 2 
dimensions, where the mechanism is simply plaquette disorder. If that were so, then 
vortices (or any other topological objects) are irrelevant at strong couplings in any 
dimension. This folklore is quite misleading, according to an argument presented in 
ref. [19], which we now review. 

Consider SU(2) lattice gauge theory at strong-coupling, and denote by U(C) the 
product of link variables around loop C. Let the minimal area of a planar loop be 
decomposed into a set of smaller areas, bounded by loops {Cj}, as shown in Fig. 
[I]. The area-law of a Wilson loop is thought to be due to "magnetic disorder", in 
which the gauge field strength fluctuates independently in sub-areas of the minimal 
surface of the loop. If this is so, then the holonomies {U(Ci)} should be (nearly) 
uncorrelated, for large areas and small (3. The test for such independent fluctuation 
in the subareas A{C{) is whether 

<nwi>=n<w.)]> (i) 

i i 

for any class function 

/••;.'/ x/a,.'/ (2) 

In fact, in D = 2 dimensions, it is easy to show that this equality is satisfied exactly. 
However, for dimensions D > 2, evaluating the left- and right-hand sides of ([!]) one 
finds for the exponential falloff on each side [19] 

e -4^(C) n l /i>>n/ie -4. m) (3) 

i 6 i 

where the inequality holds for perimeters V(C) <C Sj'P(Cj). The conclusion is that 
the holonomies U{Ci) do not fluctuate independently, even at strong-coupling, for 
D > 2. 
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If the loop holonomies in sub-areas of the loop are correlated, then where is the 
magnetic disorder required to give an area-law falloff for Wilson loops? The question 
is resolved by extracting a center element from the holonomies 



z[U(a)} = signTr[C/(Q)] e Z 2 
and asking if these center elements fluctuate independently; i.e 

<n*[tf(coi >=u<4u(c t )]> 



(4) 



(5) 



In fact, it is easy to show that they do: 
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-cjA(C) 



-oA{d) 



47T 



(6) 



Thus, magnetic disorder is center disorder in D > 2 dimensions, at least at strong 
couplings. Confining configurations must disorder the center elements z, but not the 
coset elements, of SU(2) holonomies U(Ci). Again, the only configurations known to 
have this property are center vortices. 

We now return to our conjecture that vortices are stable saddlepoints of a long- 
range effective action. Actually, there are various ways of integrating out the smaller- 
scale fluctuations, to obtain an effective action at a larger scale. One simple approach 
is to superimpose, on a lattice of spacing a with link variables denoted U, a lattice of 
spacing La with links denoted V. An effective action for the lattice with the larger 
spacing can then be obtained from 



cxp 



S eff [V]} = [DUH5[VXUU..U) V -I 



S W [U] 



(7) 



where (UU..U)i> is the product of ?7-link variables along the link I' of the V-lattice, 
and Sw is the Wilson action. Obviously, all observables computed on the V-lattice 
with S e ff will agree with the corresponding quantity computed on the U-lattice using 
Sw- 

It is trivial to compute S e ff in D = 2 dimensions, and the result is 



exp[S eff [V}} 

«nE(2j+i)(/ 2 , + i(/3)) L2 x,[v(p')] 



P , j 



cxp 



cxp 



5>g 1+ £ (2j + l) 



\ 



7" = i 1 ^ 



xAv(p')} 



+ const. 



/ 



L 2 



J2xi/2[V(P')} + const. 



(8) 
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Figure 2: Arrangement of U-plaquettes in a tube, around a rectangular contour C on the 
V-lattice. 

where V(P') is the product of V-links around the plaquette P' . One might imagine 
that the action (^j) is also a good approximation in D > 2 dimensions, at least at 
strong couplings, since this action gives the correct string tension for fundamental 
representation Wilson loops. But a quick calculation of higher representation loops 
shows that @ cannot even be approximately correct for large L. A loop in the adjoint 
representation, for example, calculated on the U-lattice with the Wilson action, is 
easily seen to have an asymptotic perimeter law falloff 

W adj (C) ~ exp[-^(C)] (9) 

where 

^ = - 41 °g(l) ( 10 ) 

is the "gluelump" mass (gluon bound to a static adjoint color charge), and V(C) 
is the loop perimeter in units of the U-lattice spacing. However, carrying out the 
same calculation with the effective action (|3|), one finds instead (with V(C) again in 
U-lattice units) the erroneous result 




A* = — 4Llog — (wrong) (11) 



with an L-dependent gluelump mass. The mismatch is not resolved by including a 
few more contours in the effective action, so long as the coupling associated with each 
contour is of order ({3/4) A , where A is the minimal area of the contour in U-lattice 
units. 

In fact, what happens in D > 2 dimensions is that the effective action contains 
Wilson loops of all sizes in j = integer representations, and these loops are only sup- 
pressed by perimeter-law coefficients. This is easily seen by bringing down a "tube" 
of plaquettes from exp(Sw) in eq. (|7|), such that the tube borders a rectangular con- 
tour C on the V-lattice, as shown in Fig. 0. Integrating over all U-links in the tube, 
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Figure 3: The degrees of freedom in Sl[V,U]. 



except those which lie on contour C, will yield contributions to S e ff such as 

i{V(C)-A) 



exp 



S eff [V}} D f DU leC II S[V$(UU..U) V - I 



x±P(C)] 




I'ec 

4(7>(C)-4) 



{Xi[V(C)\ + const.) 



(12) 



This shows that ^//[V] contains adjoint (and, in general, integer) representation 
loops with perimeter-falloff coefficients. Such non-local terms introduce non-local 
correlations among SU(2)/Z 2 coset elements in loop holonomies {U(Ci)}. Truncating 
Sgf/IV] by removing these large loops will yield erroneous results for any Wilson loop 
in representation j > |. 

Our aim is to modify the prescription (|7|) for the effective action, in such a way 
that at least the leading contribution to any Wilson loop on the V-lattice is obtained 
from a local action. The strategy for doing this is to prevent the formation of closed 
tube diagrams, of the form shown in Fig. |2], bordering contours on the V-lattice. 
This can be accomplished by not integrating, in eq. (]?]), over the U-links in a cube 
of volume 2 D around each site on the V-lattice. 2 U-links belonging to these 2-cubes 
will be denoted XJ\. In order to ease the task of illustration, we will work in D = 3 
dimensions, although the extension to higher dimensions should be straightforward. 
We then have 



DV 



DV 



IJ dUt exp[S L [V,U) 

l£2— cubes 



n 

iS2— cubes 



IJ dU v ,\{5\v^UU..U) 

l" £2— cubes 



(13) 



I' 



2 Non-local terms in the effective action will still arise from closed tubes which go around the 2- 
cubes. These, however, are associated with sub-leading contributions to color screening; the leading 
contributions in 0, arising from diagrams like Fig. will now be generated by local terms. 



6 



3 


1 


4 


2 



Figure 4: Plaquette numbering convention. 




Figure 5: Three plaquette variables on the 2-cube. 

where the long-range action Sl[V, U] depends on the V^-link variables, and on the 
[/-links in 2-cubes around sites of the V- lattice, as shown in Fig. |3|. 

Now introduce, in each 2-cube, a set of plaquette variables {hij,gij} which are 
Wilson lines beginning and ending at the center of the 2-cube, and running around 
one of the plaquettes in the cube. The hij lines run around plaquettes on the faces 
of the 2-cube, and the gr^ lines run around plaquettes on the interior of the 2-cubes. 
To fix notation: The orientation of the lines is taken to be counterclockwise in 
either the iq/-plane (i = 1), the yz-p\a.iae (i = 2), or the zrc-plane (i = 3). The 
second index (J — 1 — 4) distinguishes between the four interior plaquettes in a given 
plane with the convention shown in Fig. f|, where (x a ,Xb) = (xy), (yz), (zx). Each 

line begins at the center of the 2-cube, runs to a center of one of the faces of 
the 2-cube, goes around one of the plaquettes on the face, and returns to the center 
of the 2-cube. The orientation around the /i-plaquettes is defined by a right-hand 
rule: the thumb points in an outward direction normal to the 2-cube. The first index 
z = l,2,3 refers to a face of the 2-cube in the xy, yz, zx-planes, one lattice spacing 
away from the center of the 2-cube, in the negative z, x, y directions, respectively. 
Index values i = 4,5,6 refer to faces in the xy, yz, zx-planes one lattice spacing away 
from the center in the positive z, x, y directions, respectively. These conventions are 
illustrated in Fig. ^. 

We then integrate over the U links which do not belong to the 2-cubes. Keep- 
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ing, for each type of contribution, only terms of leading order in /3, the result is 
approximately 



J DVDU exp 

4(1,-2) 



P 



J2(Tr[h]+Tr[g]) 



+2 



73 



L 2 -4 



(14) 



Coefficients fj kl = 1 if plaquette variables /ijj and hui on nearest-neighbor 2-cubes 
can be joined by a cylinder of plaquettes, bordering link I', of length L — 2 U-lattice 



spacings. Otherwise, / z ; 



ijkl 



0. 



The next step is to change integration variables from links U to plaquettes g, h. 
This change of variables on the lattice was worked out many years ago by Batrouni 
[20], and the result is simply to introduce a Bianchi constraint into the integration 



measure 



Z « J DVDhDg JJ JJ 5[Bianchi(c(ff))] 



exp 



+2 



4(L-2) 



L 2 -4 



ETr[V^VVyy 



(15) 



Each 2-cube K contains eight unit sub-cubes; the index c in eq. ( |15D labels these 
subcubes. The 5-function constraints force a certain product of three g and three 
variables on each subcube to be the unit matrix. For example, the Bianchi constraint 
for the unit sub-cube containing hu is 



Bianchi = h^h^uh^g^ 
Now expand the 5-functions in group characters 







5 [Bianchi] = ^ (2j + l)xj[hghghg] 
i=o,§,i,. 



(16) 



(17) 



3 Thcre are 36 g, h plaquette variables on the 2-cube, and 8 Bianchi constraints, leaving 28 
independent group-valued variables. Similarly, up to 26 out of 54 link variables on the 2-cube can 
be gauge fixed to the identity, again leaving 28 independent group-valued variables. 
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and integrate out the ^-variables in eq. 
both h and /3, the result is 



). Displaying only terms of low order in 



fDVDh I] U + *{j)Y.X±[{hhh) c 

J 2-cubes K { V * / c€K 



l T l E X^hhh^ihhh), 

^ / adjacent 
CIC2&K 



x exp 



4(L-2) 



Y.fi^[h\ j v v hl l v}\ 



DVD/i exp[5 L [y,/i]] (18) 

At this stage, the action ^[V, /i] resembles an adjoint-Higgs Lagrangian, albeit 
of an unconventional form: There is an SU(2) gauge field coupled to 24 unitary 
matrix-valued "matter" fields transforming in the adjoint representation. These 
matter fields, in turn, can be subdivided into gauge-singlet fields h^^, and unit- 
modulus triplet fields eV,-, where 



hij — hijfil + i 



h 2 p- 



a 



(19) 



and e • e = 1. The unimodular degrees of freedom play a role analogous to 
Higgs fields. We know from the Elitzur theorem that the expectation values of these 
fields must vanish in the absence of gauge fixing, and cannot be viewed as order 
parameters. Since the coupling of the matter fields to the gauge field is very 
weak at large L, as compared to the self-couplings of the fields to each other on 
the 2-cubes, their expectation values depend primarily on these self-couplings, and 
on the complete removal of gauge redundancy in the e-fields through the choice of a 
unitary gauge. 

We fix to a maximal unitary gauge by first transforming one of the 24 unimodular 
"Higgs" fields eV,- on each 2-cube, denoted eU, to point in the (color) 3-direction, i.e. 



eAi = e A 2 







eA3 



(20) 



This leaves a remnant U(l) symmetry, but (|20"D is not yet a maximal unitary gauge 
fixing. We then pick one other Higgs variable on each 2-cube, denoted eg, and use 
the remaining gauge freedom to fix 



&B2 = , e B \ = sm(6 B ) > 
leaving a remnant Z2 symmetry. 



(21) 
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In the unitary gauge pO| - f21"|j , the functional integral becomes 
Z ug = I n^(n)II^( n ) A (^, h B )5(e A1 )5(e A2 )5(e B2 ) exp\S L [V, h}} (22) 



where 



A- 1 ^,^) = / d 3 a 5{e a A1 )5{e a A2 )5{em) 

d 3 a S(e ljk a j 5 3k )5(e2j k a j S 3k )5(e 2jk a j (e m S lk + e B3 S 3k )) 



From the measure 



we find 



J dh = J dh d 3 e yjl - h 2 5{e 2 - I) 



(23) 
(24) 



dh A dh B A(h A , h B )5(e A1 )5(e A2 )5(e B2 ) = J dh A0 dh B0 d9 B ^l - h 2 A0 \Jl- h 2 m sin0 B 

(25) 

Let us take, e.g., h A = hu, h B = /i 44 . Then 



Y[dV(n) / d/i 1 i )0 d/i44 ) o^44v/l — ^ii,o\A — h 2 4 sin0 44 

n J 

Yl dhij(n)exp S L [V,h] 

»i^(ll),(44) 



(26) 



The final step is to integrate over the remaining h degrees of freedom in this 
maximal unitary gauge. Defining /i-expectation values 

(F[h])h = 



— ( dhu >0 dhu,od8u)fl - ^ii,o\A - hj A sin 6> 44 J[ dh {j 



n 



^(11), (44) 



1 + 2(f) £Xi[(^)c] + 2 (f) E Xi[(hhh) cl (hhh) C2 )+... 

2 — Cubes K \ \ / C^K \ / adjacent 



n 



x exp 



F[h] 



(27) 



and 



Seff[V,h] = 2(^ 



4(L-2) 



Y,fi? kl nh\ j v v hi l vi) 



^TrfvwV*] 



(28) 
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we have 



Z ug = Z h J DV (exp S ef f[V,h] 
= Z h J DVexp[S eff [V] 



(29) 



where 



Se//IX] = S e ff[V, (h)h\ + higher-order contributions 



(30) 



The higher-order contributions consist of next-nearest neighbor (and more distant) 
couplings between (h)h terms on different 2-cubes, as well as closed loops in j — 1 
and higher representations. The large closed loops again introduce certain non-local 
correlations among the SU(2)/Z 2 elements of loop holonomies {[/(Cj)}. But S^z/IX] 
also contains local one-link terms, which provide by far the largest contribution to 
the ^-invariant part of the action. 

For the purpose of determining saddlepoint configurations of S'e/zlX] we may 
neglect the higher-order, non-local terms in the action, so that 



S eff [V] aS link [V,(h) h ] + S plaq [V] 

//?\ 4(L ~ 2) 



+2(f) i:nvvvw] 



(31) 



and, for this particular gauge choice, we find 



(hu) h 
(hn)h 
\hn)h 

(^43) h 

(h 2 i)h 
(h 2 3)h 

(^5l)/i 
(^53/ h 

(h>3i)h 
{h 3 3)h 
{h%i)h 



¥ + 



= P ~ » 



!> 

v 



!> 



3?r 



0-3 



(!) 
(!) 
'(f) 



3^3 

3<Ti 

2^ 
5 <Ti 



3tt 



0-3 



'(f) 

'(f): 
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3tt 



0"3 



• f 



'(f) 
'(f)' 



2^ 
3CT1 

2^ 

3<Ti 

^3 

37T J 



3tt 



0-3 



(hi 2 )h 
(hu)h 

\h>42)h 
(hu)h 

(h 22 ) h 
(h 2 ±)h 

\hb2)h 

(hz 2 ) h 

(^34) h 
(^64/h 



v- 

V- 
V- 
V+ 
V- 
V- 
V- 



V 



'(!) 
'(!) 
'(f) 

i|(Ti 

'(f) 



'(f) 
'(f) 



3<Ti 



2^ 
3<Ti 

2^ 
3<7i 

2^ 



3tt 



0"3 



• f 



2^ 
qCTi 



3tt 



0"3 



'(f) 
'(f) 



2 _ 



371 



0-3 



(32) 
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Inserting fl32|) into Su n k we have, to leading order in (3, 



Si ink [V,(h) h }=2 £ x£ 



4 / 

(^) 2 (f )° Tr [^Hcr3K t H] + (|) 2 (^k!l4(n)^(n)] 
Tr^n^H] - Tr^n^l^n)] 



3tt7 \4/ " V37 V4 



2 12 



37ry ,fj Tr[a s V,(n)^(n)] - (J) (|j Tr^KH^W] 

+ const. | (33) 

Each term in Su n k is proportional to a component of a IMink variable in the ad- 
joint representation, and is insensitive to the center degrees of freedom. The spatial 
asymmetry of Su n k in (0) is, of course, due to the particular unitary gauge choice. 4 
We now look for saddlepoints of S e ff. It can be seen from inspection of (|3~3"D that 
Sunk is maximized by 

V x (n) = ia 2 Z x (n) , V y (n) = ia 2 Z y (n) , V z (n) = Z z (n) (34) 

where the Z^in) = ±1 are center elements. The S p i aq term is also maximized if the 
Z^ link variables are gauge-equivalent to the identity under the remnant Z 2 gauge 
symmetry. With this choice the effective action S e ff ~ Su n k + S p i aq is maximized, 
and the configuration ([34|) is the ground state, gauge equivalent to the identity. The 
fact that this ground state is unique, up to a Z 2 gauge transformation, is again due 
to the maximal unitary gauge fixing. 
Now consider the configuration 

T/ \ ~icr 2 nx>2, n 2 = l 

VJn) = < . , 

[ +ia 2 otherwise 

V x (n) = io 2 

V z (n) = I (35) 

This configuration is a center vortex, one lattice spacing thick, running in the z- 
direction. It is not hard to see that this configuration, like the trivial ground state, 
is also a saddlepoint of S e ff. In the first place, fl3"5| ) is a global maximum of Su n k, 
since the thin vortex configuration (|3~5"D differs from link variables in the ground state 



4 Despite this asymmetry, the expectation value of any Wilson loop on the V-lattice, evaluated 
in the full (gauge-dependent) effective action S e ff[V] defined in eq. (|29|), is necessarily independent 
of the gauge choice. This should be clear from the construction, where a gauge-invariant action 
Sl[V, h] is gauge-fixed, followed by integration over the remaining h degrees of freedom. 
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only by center elements, to which Sunk is insensitive. Secondly, this configuration is 
also a stationary point of the plaquette action S p i aq [21]. This is because a plaquette 
at one of its extremal values |Tr[WV^Vt] = ±1 varies at most quadratically, and 
is therefore stationary, with respect to a small variation SV of any link respecting 
the unitarity constraint (V + 5V)(V + 5VY = I. All of the plaquettes formed from 
([35]) are extremal. So the center vortex ( p5| ) is certainly a stationary configuration 
of S e ff, the remaining question is whether it is stable; i.e. whether the vortex is a 
local maximum of S e ff, in which case it is a stable saddlepoint. 
The stability issue is settled by looking at the eigenvalues of 

5 2 S e ff _ S 2 Su n k 5 2 S p l aq 

SV^nJSV^ni) ~ 5V^ ni )5V v {n 2 ) SV^SV^) 
where, from the coefficients shown in eqs. (|3T| ) and fl33]), we see that 

5 2 S lmk (pY L ~ 2)+12 S 2 S plaq f(3 

SV^nJSVvfa) V 4 / ' ^(m)5K(n 2 ) ^ 4 

The crucial observation is that for /3/4< 1 and 

4(L - 2) + 12 < L 2 (38) 

the contribution of S p i aq to the stability matrix (and therefore to the eigenvalues of 
the stability matrix) is negligible compared to the contribution of Su n k, which has 
only stable modes. Thus from the fact that the thin center vortex ([35]) is both a 
stationary point of S e ff, and a stable saddlepoint of Su n k, we can conclude that the 
vortex is also a stable saddlepoint of the full effective action S e ff when condition 
( |38l) is satisfied. 

Condition (|38|) is satisfied for L > 5 lattice spacings. It is probably no coincidence 
that this is also where the adjoint string breaks in strong-coupling lattice gauge 
theory (as can be easily verified from looking at correlations of Polyakov lines in 
the adjoint representation). It has been known for a long time that, at intermediate 
distance scales and weak couplings, the static quark-antiquark potential is roughly 
(and maybe even accurately [22]) proportional to the quadratic Casimir of quark color 
representation. In ref. [23] this phenomenon was dubbed "Casimir scaling," and the 
problems it poses for monopole and vortex theories was discussed. In ref. [16] we 
have argued that the problems with respect to the vortex theory can in principle be 
resolved by taking into account the finite thickness of the vortex, which should be 
comparable to the adjoint string-breaking length. A result of the analysis carried 
out above is that there are stable center vortices, one lattice spacing thick on the 
V-lattice, corresponding to L > 5 on the original U-lattice. This gives an estimate 
for the vortex thickness of L = 4 lattice spacings in U-lattice units. This vortex 
thickness happens to be exactly the length where adjoint string-breaking occurs in 
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the strong-coupling Wilson action, and is therefore consistent with the reasoning in 



By inspection of S e jf, we see that the action of a center vortex configuration in 
D=3 dimensions (ignoring any further quantum corrections) is 



on the V-lattice. On the other hand, the entropy of a linelike object (essentially the 
entropy of a random walk) is a constant of 0(1) times the line length. Thus, center 
vortices are stable saddlepoints of the long range effective action whose entropy per 
unit length greatly exceeds their action per unit length. These objects therefore 
percolate throughout the lattice volume, and confine color charge in any half-integer 
group representation. 

This picture has been derived in a particular unitary gauge. A natural question 
is whether the saddlepoints of the effective action would be qualitatively different 
had we chosen to gauge- fix, instead of e*n and e 44 , some other plaquette variables 
(or combination of plaquette variables) on the 2-cubes. Although an analysis of all 
possible unitary gauge choices is beyond us at present, it is easy to see that center 
vortices must be stable saddlepoints in a very large class of gauges. We first note that, 
by definition, any maximal unitary gauge must completely determine the minimal 
action configuration of the fields, up to residual Z 2 gauge transformations. Then 
a sufficient condition for center vortex stability is simply that the classical ground 
state of ^//[V] has the form of a pure gauge V^(x) = g(x)g _1 (£ + /t), and that this 
is also a maximum of the center-invariant Su n k part of the effective action. In that 
case, the effective action must have stable thin vortex solutions at large L. This 
is because the stable fluctuation modes around a thin vortex, associated with the 
center-insensitive S^nfc term, will overwhelm (at sufficiently large L) any unstable 
modes associated with S p i aq . The vortex action at a given L will always have the 
value shown in eq. (|39D above, so the entropy of the configuration will exceed the 
action at strong couplings. 

Our findings for the strong-coupling theory do not, however, prove that vortices 
also dominate the vacuum at weaker couplings; the strong and weak coupling regimes 
are separated by a roughening phase transition, and this transition prevents a simple 
extrapolation from one regime to the other. The result is, nonetheless, significant for 
continuum physics in two ways: First, it supports the very general argument that if 
the asymptotic quark potential is sensitive only to iV-ality, then the confining field 
configurations must be center vortices. Secondly, it provides an explicit illustration 
of how center vortices are stabilized, via color-screening (center-invariant) terms in 
the long-range effective action. 

Although strong-coupling methods only apply at strong couplings, the general 
approach we have advocated here should extend, at least in principle, to weak- 



ref. [16]. 




x vortex length 



(39) 
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coupling lattice gauge theory. The central idea is that if we want to extract a local 
long-range effective action from the Wilson action, then it is necessary to include 
composite operators, transforming like adjoint matter fields, in the derivation. With 
this motivation we define gluelump operators 

G m [x;U\ = Y,"m{C x ) J] U x (40) 

which, coupled to a static adjoint source at site x, create gauge-invariant eigenstates 
of the appropriate transfer matrix. The C x are paths on the lattice beginning and 
ending at x. 5 The index M specifies the time (i.e. worldline) direction of the static 
source, and any other (e.g. spin) degeneracies. The transformation from a pure gauge 
theory on a fine lattice, to a theory of gauge fields coupled to gluelump fields H M 
on a coarse lattice, could then be accomplished as follows: 



cxp 



S eff [V,H}] = f DU\{5(y v -Q V [U\) II S(H M (x') - G M [x';U])e s ^ (41) 

J V M,x> 



where x',l' denote sites and links on the coarse lattice. The expression Qv\U] rep- 
resents a suitable "fat link" function, i.e. a superposition of Wilson lines on the fine 
lattice which run between sites bounding link V on the coarse lattice. The con- 
straints imposed by the delta-functions can be softened by replacing delta-functions 
with exponentials, as in the "perfect action" approach [24]. The end result of this 
procedure will be an effective long-range action consisting of gauge fields coupled to a 
set of adjoint Higgs-like fields. Possibly this scheme can be implemented numerically 
at moderately weak couplings, along the lines of the Monte Carlo renormalization- 
group. 

In the strong-coupling analysis carried out in this article, we have seen how 
color-screening, center-invariant "Higgs" terms predominate in the long-range ef- 
fective action beyond the adjoint string-breaking scale, and stabilize center vortex 
configurations. The entropy of these configurations exceeds the cost in action, and 
vortex configurations percolate throughout the lattice. We think it likely that these 
important features are not specific to strong couplings, and also characterize the 
effective action of lattice QCD in the continuum limit. 
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